Where does a cohesive granular heap break? 
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In this paper, we consider the effect of cohesion on the stability of a granular heap and compute 
the maximum angle of stability of the heap as a function of the cohesion. We show that the 
stability is strongly affected by the dependence of the cohesion on the local pressure. In particular, 
this dependence is found to determine the localization of the failure plane. While for a constant 
cohesive force, slip occurs deep inside the heap, surface failure is obtained for a linear variation of 
the cohesion on the normal stress. Such a transition allows to interpret some recent experimental 
results on cohesive materials. 

PACS numbers: 61.43.Gt, 61.43.Gt, 45.70.Cc 
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I. INTRODUCTION 

The economic impact of particle processing is enor- 
mous. Better methods for the design and synthesis of 
unit operations involving divided solids have been iden- 
tified as a critical need, especially for pharmaceutical, 
agrochemicals and specialty chemicals. One of the im- 
portant features in industrial processes is to better un- 
derstand the transition between a static granular medium 
and the avalanche process. Avalanches in non-cohesive 
granular media have been extensively investigated 0] • 
A common characteristic of these studies is that granular 
motion occurs in a relative thin boundary layer (around 
ten grains) at the surface independent of the size of 
the sample. On the other hand, recent experiments have 
explored the relatively new subject of cohesive granular 
media such as "humid granular" media. The presence of 
capillary bridges between the grains generate adhesive 
forces which str ongl y affect the stability of the samples 
H H S H H Efl 113 ■ These capillary bridges either 
originate in small amounts of added fluid, and this sit- 
uations corresponds to "lightly wet granular media" or 
either are created by a condensable vapor in the atmo- 
sphere HEI!. In the latter case the vapor is in chemical 
equilibrium with the interstitial liquid bridges. These ex- 
periments concern the so-called "moist granular media" . 

As shown by all these experiments, the behavior of 
humid granular media strongly departs from the dry ma- 
terials and many new features appear. In a series of ex- 
periments on "lightly wet" granular media, Tegzes et al. 
|12j . measured the angle of repose of a granular media 
(by the draining crater method) as a function of the liq- 
uid content and of the size of the system. They observed 
an increase of the angle of repose with the liquid content. 



More precisely, they obtained three regimes: i) at low liq- 
uid content, the "granular regime", the avalanche occurs 
at the surface and the angle of repose does not depend on 
the size of the system; ii) at intermediate liquid content, 
the "correlated regime" , the avalanche takes place more 
deeply in the heap, the angle of repose depends on the 
size of the system; iii) at high liquid content, a kind of 
"plastic flow" is observed. This various regimes are not 
accounted for by previous theoretical analysis. 

In numerical and analytical modelizations, cohesion is 
usually taken into account as a constant adhesion force. 
As an example, Forsyth et al. ^3 and Valverde et al. ^3 
used a constant van der Waals force to interpret their ex- 
perimental results, Olivi-Tran et al. EH used a constant 
capillary force or a constant force due to a solid bridge 
in particle dynamics (PD) simulations, Albert et al. 0] 
used a constant capillary force in a geometrical model 
and Nase et al. [13, Mikami et al. [Ig used a constant 
capillary force in PD simulations. 

In this paper we show that the assumption of a con- 
stant cohesion is not justified in many practical cases. 
We then show that the dependence of the cohesion on 
normal stress, usually omitted in the literature, is a key 
point in determining the stability of a cohesive granular 
heap. This effect leads to a dependence of the stabil- 
ity on the system size, as found experimentally. To this 
end, we shall use a continuum analysis [l7| to study the 
stability of cohesive sandpiles taking into account a non 
constant cohesive stress in the pile. The paper is divided 
in two parts. We briefly first show that in a many sit- 
uations, the adhesion force between the grains depends 
on the normal stress. Then, in the second part, we show 
how this relationship affects the stability of a heap and 
the localization of the slip plane upon failure. 
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II. CONTACT FORCES BETWEEN GRAINS 

The adhesion force between two ideal spheres due to 
solid-solid interactions has been calculated in the case of 
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an elastic contact a long time ago 0, 0, |2(j : 

F s _ s = /7r 7 i? (1) 

where R is the sphere radius, 7 is the solid/interstitial 
medium surface tension (7 = "/sg if the spheres are in a 
gas atmosphere, and 7 = 7 sl if they are immersed in a 
liquid), and / is a numerical factor between 1.5 and 2. 
The lower value of / = 1.5 corresponds to the Johnson- 
Kendall-Roberts calculation (JKR) valid when the at- 
tractive forces are strong enough to deform the spheres 
surfaces. The higher value / = 2 corresponds to the 
Derjaguin-Muller-Toporov (DMT) calculation, which is 
valid in the opposite limit of weak attraction and rigid 
solids. The main hypothesis in obtaining Eq. ^is that 
surfaces are ideal. An obvious feature of this result is 
that adhesion does not depend on the pressure applied 
on the surfaces before to separate them [2l|. However 
measured values of the adhesion force generally departs 
from the predicted one [2^. As we briefly discuss in the 
next section, the effect of surface roughness leads to qual- 
itatively different results, involving a dependence of the 
adhesion force on the normal load. 



A. Pressure depending adhesion between two 
grains 

Restagno et al. (2^ have recently measured the adhe- 
sion force between two weakly rough surfaces (a sphere 
and a plane) using a surface forces apparatus. The "pull- 
off" force -Fkdh to separate two surfaces has been found to 
depend on the a loading normal force, F, with a scaling: 
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(a) (b) 

FIG. 1: (a): Contact between two perfect spheres of radius 
R, Young modulus E. 

(b): Contact between two rough particules of radius R. 

The key-point to interpret this result is the roughness 
of the surfaces as pointed out 50 years ago by MacFarlane 
and Tabor |23j: the surfaces are in contact only at the 
tip of their asperities, and the real contact area ^4 r is 
much smaller than the apparent contact area ^4 app (Fig. 
n To obtain the adhesion force they used Eq. ^ with 



an effective surface tension 7e g = 7^4 r /^4 app . On one 
hand, different ways can be used to calculate A r . The 
first way is to consider that the pressure in the molecular 
contact area is so high than it can reach the yield stress 
of the material 23] . An other way has been proposed by 
Greenwood [24| considering the statistical distribution of 
heights of the surface asperities. This two models lead 
to a simple proportionality between the real contact area 
A r and the normal load: 

A r ~ F (2) 

On the other hand, the calculation of the apparent con- 
tact is the elastic contact area calculated by Hertz [25| : 

A ap p ~ (FR/E*) 2/3 , with E* = E/[2(l - v 2 )] the re- 
duced Young modulus of the material, defined in terms 
of the Young modulus E and the Poisson modulus v of 
the material. Gathering results, one obtains eventually 
7 c ff ~ F 1 / 3 , so that Eq. ^is compatible with the observed 
experimental scaling. 

In the previous experiments, the cohesion between two 
rough surfaces due to the direct interaction between the 
surfaces has been considered. However, adhesion might 
also originate in the presence of a cohesive "binder" in 
between the surfaces. In many practical situations, the 
presence of liquid bridges connecting the grains leads to 
strong adhesive forces between the grains 0, IE IE IE 
9]. When the material is placed in an atmosphere of 
a condensable vapor 0, EJ , the size of the liquid bridge 
is fixed via chemical equilibrium by the value of external 
humidity. The size of such bridges is very small, around a 
few nanometers, and one expects that the liquid bridges 
will condense only in the regions close to real contact: 
One expects therefore that F co h ~ A T , the area of real 
contact area (see e.g. Q for a more complete discussion). 
Since as discussed above in Eq. [3 A r is proportional to 
the external load, one gets eventually: F a dh ~ F. This 
picture is however valid when the size r of the bridges 
is lower than a typical scale of roughness, £r: In the 
opposite case, r > £r, roughness is not pertinent any 
more and the capillary force reduces to a contant value 
2ir-fR, with R the radius of the beads, and 7 the liquid 
vapor surface tension. 

To sum up, the adhesion force between two surfaces 
can be written: 

Fadh ~ F n (3) 

with various expected exponents n. For humid granular 
materials, a transition from from n = 1 to n = is 
expected, depending on the relative size of the capillary 
bridges (i.e. on the added amount of liquid) compared 
to roughness : For small amount of liquid, the value n = 1 
is expected, while large amounts of liquid correspond to 
n = 0. On the other hand, for dry or liquid saturated 
materials in which cohesion is expected to originate in 
direct (e.g. van der Waals) interactions, a value n = 1/3 
is predicted. 

We now derive the local cohesion stress from the ad- 
hesion force. 



3 



B. from adhesion force to local cohesion stress 



III. FAILURE OF A COHESIVE HEAP 



The problem of the derivation of the interparticle 
forces F from the bulk stress a is an old problem. Quin- 
tanilla et a/.|lll| have shown, comparing direct adhesion 
between particules using an AFM and cohesion stress 
measured by determining the tensile stress of a mate- 
rial, that adhesion increases with the normal load in the 
material and that a correlation between the microscopic 
contact force and a macroscopic stress measurement can 
be performed with a good precision. 

For a system of hard monodisperse spherical particles 
with a random isotropic packing, a continuum-like pic- 
ture leads to a simple linear relationship between force 
F and stress a Hal : 



<f>k 
4vri? 2 



F 



(4) 



where R is the radius of the particles, k is the coordina- 
tion number, which is defined as the average number of 
contacts per particle, and is the volume fraction. Such 
a relationship has been widely used in the literature for 
powders or wet granular media [Til [271 • From Eq. 0] we 
can deduce the relationship between the adhesion force 
-Fadh and the adhesive stress c ac jh in the granular medium: 



Cadh 



4irR 2 



Fadb 



F 
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(5) 



With the same arguments, we can deduce the relationship 
between the normal load and the normal stress: 



4nR 2 ~ R 2 



(6) 



Assuming a general relationship of the form of Eq. 
one gets therefore an adhesive stress which depends on 
the normal stress as shown in Fig. |21 



Cadh — Cq 



00 



(7) 



where n, <7o and cq characterize the properties of the 
material. This relationship is our general starting point 
to study the stability of a cohesive heap. 



n > I 

"Apparent non-cohesive regime" 



n = l 

"Granular regime" 




0<n< 1 

"Size dependent avalanche regime" 
» = 

"Classic adhesion regime" 



FIG. 2: different regimes of cohesion in a granular heap. 



The most simple approach in the continuum limit to 
study the failure of a granular material is the plastic the- 
ory of Mohr-Coulomb. The validity of the continuum 
approach has been widely discussed. Nevertheless con- 
tinuum approaches allow to predict the Green's function 
of a granular layer [3l| and is widely used in soil me- 
chanics [rH W\ . We consider the stability of a cohesive 
granular heap, characterized by an adhesive stress de- 
pending on the normal stress, as given in Eq. [3 in which 
we suppose the Mohr-Coulomb criterion valid. Our aim 
is to locate in such a material the " slip plane" , where 
failure occurs. The basis of our analysis is the Coulomb 
criterion: a granular material is stable if for each surface 
inside the material the following inequality is obeyed: 



t < /ACT 



(8) 



where r is the shear stress and a is the normal stress. For 
a cohesive material, we assume that this condition can be 
generalized by adding in Eq. [SJ the adhesive stress c a( un 
as defined in the previous section, to the normal stress a 



We consider the geometry depicted in Fig. [3] A gran- 
ular heap, with height H, makes an angle 9 with the 
horizontal. The heap is supposed to be invariant in the 
direction perpendicular to the figure. We follow and gen- 
eralize the approach given in Ref. [l^ to the case of a 
cohesive material with the law of cohesion given in Eq. \7\ 
In the simplified description of Ref. ^3] j the sn P surface 
is assumed to be planar, making an angle a against the 
horizontal. 




FIG. 3: Geometry of the granular heap. 6 is the actual slope 
of the heap, while a locates the position of a possible failure. 



The force balance in the direction perpendicular and 
parallel to the slip plane leads to the following conditions: 



P sin a N = P cos ot + F t 



adh 



(9) 



where P is the weight of the granular material above the 
slip plane; c = F/S s n p and a + c a dh = N/S s y lp are the 
normal and tangential stress along the possible slip plane 
with area S s y lp . In these variables, the stability criterion, 
Eq. [S] writes F < fiN. A simple geometric calculation 
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leads to the relationships [l7| : 



P 



S s Up — 



/ 1 



2 \ tan a tan( 



sma 



(10) 



with L the length of the heap in the invariant direction 
(perpendicular to the figure in Fig. Gathering these 
results, one gets the following stability criterion: 



( 1 



\tana tan 



1 \ , . n . c a dh H 

(sin a — n cos a) < \i — - 



sma 



(11) 

Introducing the angle 9q defined as tan 9q = M> this rela- 
tion can be adequately rewritten as: 



sm(8 — a) sm(a — 9q) < 2sin6'o sini 



, C a dh 

pgH 



(12) 



Now, one has to introduce the dependence of the co- 
hesive stress c a dh as a function of the normal stress, as 
discussed in the previous section. We shall use the gen- 
eral expression given by Eq. [3 Using a — P cos a/ S s ii P , 
one gets after some algebra: 



Cadh — Cq 



/ pgH\ n / sin(9 — a) cos a 



smt 



V 2a / 

Introducing the cohesion parameter £ co h defined as: 

> cq ( pgH \ n 

Ccoh ~ pgH\2aJ ' 

and the function f[a] defined as: 



„ r n (sin(0 - a)) 1 "" 

/[<*] = L--i — sin (« - d °) (1 5 ) 



(13) 



(14) 



(cos a)™ 



Eq. H2l rewrites: 



f[a] < Ccohsm6>o(sin0) 



n-l 



(16) 



The heap is therefore stable at an angle 9 if this inequality 
is verified for all possible a (with < a < 9). In the 
opposite case when this inequality is not verified for some 
values of a, the slip plane corresponds to the value of a 
for which this inequality is "first" violated. Note that for 
a non-cohesive material, C co h = 0, the maximum angle of 
stability is simply (9o- 



A. n < 1 case 

First if 9 < 9q, the function f[a] is always negative and 
the stability condition is trivially verified: For 9 < 9q the 
heap is always stable. 

The opposite case 9 > 9 is more complex. A typi- 
cal plot for this function is then given in figure 0] This 
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FIG. 4: Plot of f[a] as a function of a, for a specific choice 
of parameters: 8 = 60°, 8o = 30° and n = 1/3. The dashed 
line locates the position of the maximum a m . 



function does exhibit a maximum for a value a m , ver- 
ifying df/da(a m ) = 0, and depending of the heap an- 
gle 9. According to criterion ^3 the heap is stable if 
/[c>! m ] < Ccohsin6'o(sin6')" _1 . The maximum angle of 
stability 9 nl thus is solution of the implicit equation: 



/W = CcohSin6'o(sin( 



(17) 



A straightforward calculation shows that a m obeys the 
following relationship: 

sin(2ct m — 6*o — 9) cosa m = ncos(2a m — 9) sin(a m — 9q) 

(18) 

It is interesting to consider the two limiting cases n = 1 
and n = 0. 

• for n = 1, the solution of Eq. 1181 is: a m = 9; and 
9 m is found to obey: 



sin(fl m - go) 
cos 9 m sin 6*o 



= c 



coh 



(19) 



This relationship can be rewritten in a more explicit 
form as : 



tanf 



Ml +Ccoh) = p ( 1 + 



(20) 



with fi — tan# - This case presents two important 
characteristics: (i) failure occurs here at the sur- 
face of the heap, as emphasized by the relationship 
a m = 9; (ii) whatever the heap height H, there is 
always an effect of cohesion on the maximum angle 
of stability. This regime has been called the "gran- 
ular regime" according to Tegzes experiments in 
which a cohesion effect on the avalanche angle is 
found without dependancy on the heap size. In 
Tegzes experiments, this regime is found at low liq- 
uid content. In sectionlul we showed that the value 
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n = 1 is indeed expected when cohesion results of 
small capillary bridges between surfaces asperities. 

• for n = ("classic adhesion regime" i.e. con- 
stant adhesion force), the solution of Eq. 1181 is: 

Q^m = ($o + $)/2, and 9 m verifies the equation: 

1 - cos(0 m - Op) _ _ c 
2 sin 9 sin 0q pgtl 

This is a classical result, as obtained e.g. in Ref. 
[l7| . It is important to note that in this case, (i) the 
heap fails deep inside the material, since a m < 9 m ; 
(ii) the maximum angle of stability 9 m depends on 
the height of the heap H. Such a dependence is in 
fact expected when one realizes that for n — 0, a 
"capillary length scale" can be defined on dimen- 
sional grounds : i cap = c / pg (see eg Eq. El ■ This 
capillary legnth gives the size of a macroscopic ele- 
ment of the granular medium whose weight equals 
the adhesion force which act on it. For instance 
it maybe the maximum size of a powder aggregat 
which can remain stuck under a horizontal surface. 
Eq. shows that the maximum angle of stability 
of the heap reduces to 9q when the size of the heap 
is very large compared to the "capillary length" . 

• For n in between these two limiting cases, one has 
to solve numerically Eqs. El and El a m lies m 
between 9 and (9 + 9q)/2. We show in Fig. a 
typical result for the dependence of 9 m and a m as 
a function of Ccoh for n = 1/3. The dependence of 
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FIG. 5: Dependence of the maximum angle of stability, # m 
(solid line) and failure plane location, a m (dotted line), as a 
function of the cohesion f co h (angles are here given in degrees). 
Physical parameters are: 9o = 30°, n — 1/3. The dashed line 
is the approximate solution, as given by Eq. 1221 

9 m in the limit of small cohesion can be computed 
analytically. For Ccoh = 0, one has a m = 9 m = 9q 
and one may expand the angles around this values 



for small Ccoh: «m = 9 a + 5a and 9 m = 9 a + 59. 
First linearizing Eq. 118(1 . one gets the relationship 
5a — 59/(2 — n). Introducing this condition into 
Eq. IT7|) . one gets eventually: 

0m - 9 = 7 (22) 

with 7 = (2 n)/(l 

n )(l-n)/(2-n) sin o ( cos o )n/2-n Inserting the 

H dependence of Ccoh, as defined in Eq. El one 
gets therefore 

9 m - 9 a oc H 5^ (23) 

This shows in particular that the dependence of 
the maximum angle on the heap height is a direct 
measure of the power n, characterizing the nor- 
mal stress dependence of the cohesion stress. This 
regime is called the "size dependent avalanche" 
regime. 

B. n > 1 case 

In this case, it is easy to show that whatever the cohe- 
sion, the maximum angle of stability is 9q. This regime is 
called the "apparent non-cohesive" regime. First as 
in the previous case, when 9 < 9q, the condition of sta- 
bility, Eq. Hfijl . is trivially obeyed since the function f[a] 
is negative. Now, for 9 > 9 , f[a] goes to infinity when 
a — » 9. The heap is therefore always unstable, whatever 
the cohesion. As a result, 9 — 9q is the maximum angle 
of stability, independently of the cohesion £ co h 

IV. DISCUSSION 

In this paper we have studied the effect of an adhesive 
stress on the localization of the slip plane in an avalanche 
process. More precisely, we have shown that in a lot of 
practical case, the adhesive stress in a material depends 
on the normal stress. We have calculated, using a Mohr- 
Coulomb analysis, the internal angle of slip a m . 

A few conclusions can be drawn from these results: 

• First, the location of the slip plane, here defined 
through a m , does also strongly depend on this func- 
tional dependence: While for a constant cohesion 
(i.e. independent of the normal stress) the heap 
slips deep inside, only surface slip is expected when 
cohesion is linearly related to normal stress. 

• Second the maximum static angle, 9 m , does depend 
on cohesion via the dimensionless parameter, Ccoh, 

defined as Ccoh = i^iof) ' ^ n i m P or tant re- 
mark is that for n / 1, the maximum static an- 
gle depends on the height of the heap, H. On the 
other hand, for the specific value n = 1 - which 
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is expected in some physical situations (see discus- 
sion in section 111 Aj l- this dependence disappears 
and 9 m is an intrinsic property of the material, in- 
dependently of the geometry. 

• Eventually, a change of regime in the cohesion, e.g. 
a change from n — 1 to n = as discussed in section 
111 Al will not only modify the dependence of 6 m on 
cohesion, but more dramatically, it will change the 
localization of the slip plane: For example, while 
for n — 1 slip occurs at the surface of the the heap 
(a m = m ), it will fracture in the interior of the 
heap for n = (a m < 9 m ). In other words, any 
change of slip behaviour reflects a transition of co- 
hesion regime. 

This features provide a framework to explain the re- 
sults of Tegzes et al.. In view of the previous results, 
it seems natural to interprete the different regimes ob- 
served in the experiment, as an indication of a change in 
the cohesion regime. At low liquid content, the avalanche 
takes place at the surface of the heap and the angle in- 
creases linearly with the adhesion. With our approach, 
this would correspond to the n — 1 regime, which is in- 
deed expected at low liquid content (where the capillary 



bridges do not fill the full interstitial space between the 
grains). At higher liquid content, a change to a regime 
with n < 1 is expected. This might explain the various 
observed experimental features, such as the size depen- 
dence. 

As a conclusion, we hope that this work will motivate 
further experimental investigation on the stability of co- 
hesive granular materials. The present results suggest 
that a careful determination of the stability properties 
and of the failure plane localization yields information 
on the cohesion forces between grains. In particular the 
heap height dependence of the maximum angle of stabil- 
ity should provide a direct measure of the cohesion prop- 
erties. Another interesting geometry is the cylindrical 
bunker ("Janssen's problem"), in which cohesion effects, 
as discussed here, should play a particularly important 
role. 
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